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Abstract 

Approximate relations among transverse momentum dependent quark distribution 
functions are established in the framework of the QCD parton model. The validity 
of such results survives QCD evolution effects, owing to the Politzer theorem on 
equations of motion. Furthermore the model fixes an energy scale, involved in the 

parametrization of the correlator, which determines the Q 2 dependence of the 
azimuthal asymmetries in inclusive reactions. Some of the present data - the cos2cj) 
asymmetry in unpolarized Drell-Yan, the sin2<fr single spin asymmetry in 
semi-inclusive deep inelastic scattering (SIDIS) and the coscj) asymmetry in 
unpolarized SIDIS - support model predictions. Further measurements of SIDIS and 
Drell-Yan asymmetries are suggested, in particular the double spin asymmetry in 
semi-inclusive deep inelastic scattering, which allows to determine approximately 

the proton transversity. 



PACS numbers: 13.85.Qk, 13.88.+e 



1 Introduction 



The quark-quark correlator [1, 2] is a fundamental tool for calculating cross sections 
in inclusive high energy reactions. It may be parametrized according to the Dirac 
algebra, using the available vectors. We distinguish between the common correlator, 
which depends only on the longitudinal fractional momentum, and the transverse mo- 
mentum dependent one. In the former case the correlator consists of 9 components [3], 
while in the latter case the number of independent functions is 32 [4, 5] - two for any 
Dirac operator - of which 23 are washed out upon integration over transverse momen- 
tum. These last functions contain the dimensionless four-vector 



where p± is the transverse four-momentum of the active quark with respect to the 
hadron and /j,q an undetermined energy scale, introduced for dimensional reasons. 
We shall fix this parameter by comparison with the limiting case of the QCD par- 
ton model, which provides also useful approximate relations between different soft 
functions. Some predictions of this model are tested against data and further mea- 
surements are suggested. 

First of all we introduce the quark-quark correlator, $, discussing its normal- 
ization. Secondly, we parametrize <£> according to the Dirac components. Then we 
write the density matrix of a quark in the approximation of the QCD parton model. 
Comparison with the general parametrization of $ allows to determine fio and implies 
approximate relations among transverse momentum distributions. Lastly we compare 
our results with present data and suggest for them further experimental tests. 

2 Parametrization of the correlator 

The correlator is involved in a typical calculation of inclusive cross section - for exam- 
ple, deep inelastic scattering (DIS), semi- inclusive DIS (SIDIS), or Drell-Yan (DY) - 



V± =P±/Ho, 
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Here dT is the phase space element, a the fine structure constant, Q 2 the virtuality of 
the exchanged photon and T the flux factor. Moreover IS" and W^ u are respectively 
the leptonic and hadronic tensor. In particular, the hadronic tensor reads, at zero 
order in QCD, 

W» u = cJ2e 2 a [ d 2 p ± Tr [<S> a A (x a , p ± )-f* b B (x b , - p ± )-f] . (3) 

Here c = 1 and the flavor b coincides with a for DIS (and SIDIS), while c = 1/3 and 
b = a for DY. Each correlator is normalized so as to reduce to the usual spin density 
matrix in absence of parton-parton interactions. For a nucleon we have 

$ = $ e + *o, (4) 

where $ e is even under time reversal and $ D is odd under the same transformation. 
$ e reads, up to and including twist 3[2, 6, 5], 

$ e ~$(2) + $(3) ; (5) 



with 



V 1 

+ \(Mi L + \ ± hiMji±^ + }} (6) 



and 



*?> = ^(/ ± + A^75 + A^75)^ + ^A ± ^7 5 [^^J 

+ -xM{e + ^75^ + -(A/i L + A ± /i T ) 75 [^_,^ + ]}. (7) 

On the other hand, the twist-2 component of <3> amounts to [7, 8] 

$ ° " ^^wWiSI + ^ i ^i.f'+]}- (8) 

In formulae (6) to (8) we have used the notations of refs.[2, 9] for the "soft" functions, 
but with different normalizations, as we shall specify below. They are functions of 
x and of p 2 ± , where x = p + /P + and p± are, respectively, the longitudinal fractional 
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momentum and the transverse four-momentum of the quark. We have chosen a 
frame where, in light cone coordinates, the four-momentum of the nucleon is P = 
(P + ,P _ ,0), the quark four-momentum is p = (p + ,p~ , p_i_) and p±_ = (0,0, pj_). n± 
are lightlike vectors, such that n + ■ n_ = 1, whose space components are directed 
along (+) or opposite to (-) the nucleon momentum. Moreover P 2 = M 2 and 

S = X^ + S ± (9) 

is the Pauli-Lubanski (PL) vector of the nucleon; one has S 2 = — 1, A = —S ■ n and 
n = 1/^/2(1, —1,0) in the nucleon rest frame. Thirdly, 

v = -L P . n _, x± = -s- v± , (io) 

r)± being given by eq. 1. Lastly, the energy scale /to, encoded in rj±, has been 
introduced in such a way that all functions involved in the parametrization of <£> have 
the dimensions of a probability density. This scale - defined for the first time in 
ref.[10], where it was denoted by - determines the normalization of the functions 
which depend on therefore /io has to be chosen in such a way that these functions 
may be interpreted just as probability densities. We can reasonably assume this 
parameter to be independent of the perturbative interactions among partons. 

We conclude this section with a comparison between the normalization of our 
functions and the one of Mulders et al.. We distinguish between the functions which 
are not multiplied by rj± and those which are multiplied by that vector. As regards 
the former category, we have, e. g., 

A = f?/x, giL = g? L /x, h 1T = h? T /x, (li) 

where the superscript m refers to the Mulders convention. Concerning the second 
group of functions, the relations are of the type, for example, 

grr = 9? T M/M, h\ L = h^M/i^x). (12) 
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3 Density matrix in QCD parton model 



Now we write the density matrix of a confined quark, but free of interactions with 
other partons. To this end, we inspire to the one of a free spin- 1/2 particle, i. e., 

Pf = l(t + m)(l + l5 n (13) 

where p, m and S' are respectively the 4-momentum, the mass and the Pauli-Lubanski 
4- vector of the fermion. We rewrite eq. (13) as 

Pf = \(j + m)(l + 75^1 +75^1), (14) 

where Sy = X'p/m, S' ± — S' — S'^, X' — p • s and p and s are unit vectors in the 
directions, respectively, of the quark momentum and of the spin in the particle rest 
system. In order to generalize eq. (14) to the case of interest, we take into account 
three elements: 

a) the internal structure of the nucleon is sensitive both to longitudinal and to 
transverse polarization, in an independent way; 

b) the spin of a massive particle has to be defined in its rest frame, therefore the 
PL vector of the quark does not coincide with the one of the nucleon[ll]; 

c) the quark can be treated as if it were on shell, owing to the equations of motion; 
we shall discuss this point below. 

Then 

Pc = \p> + m q ) [q(x,p 2 ± ) + Aq(x,p]_)-f 5 ^ + A T q(x,p 2 j_)-f 5 S q ± ] . (15) 

Here and Sj_ are defined in such a way that, in the quark rest frame, they coincide 
respectively with S\\ — S — S± and S± in the nucleon rest frame. Therefore 

SI = X(JL - r} ± ) + 0(rjl), Sl = S ± + X ± ^- + 0(rj 2 ± ), (16) 

Ulq 17lq 

with rj ± = p±/V and Aj_ = —S-rj ± . Moreover, we have set 

q(x,p 2 ± ) = q+(x,p 2 ± ) +q-(x,p 2 ± ) = q^(x,p 2 ± ) +qi(x,p 2 ± ), (17) 
Aq(x,p 2 ± ) = q + (x,p 2 ± ) - q-(x,p 2 ± ), (18) 
A T q(x,p ± ) = qi(x,p 2 ± ) - qi(x,p 2 ± ). (19) 



Here q±(x,p 2 ± ) and qm)(x,p 2 ±) are the probability densities of finding a quark with, 
respectively, a positive (+) or negative (-) helicity and a positive (f) or negative (1) 
transversity, the sign of the latter being determined by the sign of the scalar product 
between the quark transversity and the nucleon transversity. Eq. (15) can be obtained 
as a limiting expression[12] of the correlator in a gauge theory, i. e., 

$ = j &(p;P,S)dp-, (20) 

where the matrix elements of $'(p; P, S) are defined as 

%(P-,P,S) = J ^e^(P,S\^(0)£(x)Mx)\P,S). (21) 
Here ip is the quark field and \P, S) is a state of the nucleon. Moreover 

C(x) = Pexp[-igA P (x)], with A v (x) = [* \ a A a (z)dz», (22) 

Jo ^ 

is the gauge link operator. Here "P" denotes the path-ordered product along the 
integration contour V, \ a and A a ^ being respectively the Gell-Mann matrices and the 
gluon fields. The link operator depends on the choice of V, which has to be fixed 
so as to make a physical sense. However, for our aims we can neglect the details 
of the contour. Indeed, in the limit for g — > 0, that is, for noninteracting quarks, 
£(x) tends to 1. Moreover, in that limit, the quark may be treated as an on-shell 
particle, as shown also by Qiu[13] via equations of motion[14]; see also ref.[15]. It can 
be shown[12] that in this limit $ tends to p c . Another important consequence of the 
equations of motion is that they survive renormalization, therefore eq. (15) is a good 
approximation to the correlator even if QCD evolution of the distribution functions 
is taken into account. 



4 Approximate equalities among "soft" functions 

Now we compare the density matrix (15) with the general parametrization (6) of the 
twist-2, T-even component of $. To this end we consider projections of both matrices 
over the various Dirac components, i. e., for a given Dirac operator T, 

$ r = iTrr$. (23) 
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First of all, T = 7 + , 757+ and 757+7; (i — 1,2) yield 

h = Q, 9il = Ag, hiT = A T q (24) 

and 

h iL ~ ~^Ag, # 1T « /i^ « ^A T g. (25) 

The last equalities hold approximately in the limit of m g = 0. In order to determine 
fj,o, we observe that the functions g±T, h\ L and hj~ T , involved in formulae (25), are 
twist 2, therefore they may be interpreted as quark densities, as well as Aq and A T q. 
For example, gyr is the helicity density of a quark in a tranversely polarized nucleon. 
Therefore it is natural to fix /i in such a way that such functions are normalized like 
Aq or A T q. This implies, neglecting the quark mass, 

^ = V = -^=p-n_. (26) 

This result differs from the treatments of previous authors[2, 16], who assume /j,q = M. 
Moreover, it is interesting to consider also the projections over twist-3 operators, in 
particular r = 7$ (i — 1, 2). This yields 

f L » fi = q, (27) 

which is known as the Cahn effect [17, 18]. We have neglected quark-gluon interac- 
tions[2], which could modify, in principle, equalities (24), (25) and (27), in particular 
the last one, concerning a twist 3 operator. However, as we shall see, comparison with 
data, wherever possible, suggests that such equalities are approximately verified. On 
the contrary, the projection over T = 757^ (i — 1,2) yields (after integration over pj_) 

9t(x) = ^(x). (28) 

In this case the contribution of the QCD parton model is very small: m q is negligi- 
ble for u- and ci-quarks, while for s-quarks h\ is predicted to be small, because sea 
quarks are produced mainly by annihilation of gluons, whose transversity is zero in 
a nucleon. Therefore the contribution of quark-gluon interactions, neglected in our 
model, becomes prevalent in this well as for T = 1 and 757+7-, corresponding 
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respectively to e and h^. Such interactions will be discussed below. Eq. (28) - simi- 
lar to some expressions given in the literature [19, 2] - establishes a relation between 
transversity and transverse spin. Indeed, the two quantities are related to each other, 
as can be seen also from eq. (14) for a free fermion. But, unlike transversity, the 
transverse spin operator is chiral even and does not commute with the free hamilto- 
nian of a quark[3]: in QCD parton model it is proportional to the quark rest mass, 
which causes chirality flip. 
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Figure 1: Azimuthal asymmetry in unpolarized DY: the asymmetry parameter v vs 
p = |q_i_|/Q- Data are taken from the first two refs.[20]: circles correspond to = 
16.2 GeV, squares to ^/s = 19.1 GeV and triangles to ^fs = 23.2 GeV. The best fit 
is made by means of formula (30), taking into account eq. (29). A = 1.117. 

Applying eq. (26) to SIDIS and DY yields 

/io = Q/2. (29) 

Moreover, results (24) to (27) can be extended, with some caution[4], to the fragmen- 
tation correlator. 

It is worth comparing our approach to previous ones. Kotzinian[10] starts from the 
approximate expression of the density matrix for a free ultrarelativistic fermion and 
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adapts it to the case of a quark in the nucleoli. He parametrizes the density matrix 
with the 6 twist-2, T-even functions that appear in the parametrization (15). Similar 
results are obtained by Ralston and Soperfl] and by Tangerman and Mulders[9]. 
The difference with our approach is that those authors do not take into account 
the Politzer theorem, which implies relations among the 6 soft funtions. In any 
case, eq. (15) turns out to give the same result as found by Ralston and Soperfl] 
as regards the parametrization of the DY hadronic tensor in the limit of parton 
model. Indeed, according to our approximation, the symmetric hadronic tensor for 
two polarized protons consists of just three coefficients, corresponding, respectively, 
to the convolutive products q <g>q, Aq ® Aq and Axq <8> Axq. 

A remark is in order. Eqs. (24), (25) and (27) have been deduced in the frame- 
work of the QCD parton model. They survive QCD evolution, but are modified by 
quark-gluon interactions. In particular, the Politzer theorem implies [15] that such 
interactions modify those approximate equalities by terms of order g 2 /Q 2 . This pre- 
diction could be tested, for example, in the case of gr, by performing measurements 
at different Q 2 values. 



5 Phenomenological analysis 

Now we test the relations just found against results from three different experiments. 

1) Unpolarized Drell-Yan with a fixed transverse momentum of the final muon 
pair. In this case the momentum of each muon presents an azimuthal asymmetry 
with respect to the plane containing the initial beam direction and q^. This az- 
imuthal asymmetry is essentially of the type l/2h / sin 2 9cos2(p[20], where 9 and are 
respectively the polar and azimuthal angle of the fi + momentum in the Collins-Soper 
frame[21], defined in the center of mass of the pair, while v is a dimensionless coef- 
ficient. This asymmetry is interpreted [20] in terms of quark-ant iquark annihilation 
into a virtual photon. In the formalism of the correlator[2, 8, 12] the asymmetry 
is described[23] as a convolutive product of the transversities h\ and hf of the two 
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active partons in the initial hadrons, providing for v an expression of the type [12] 

v = Aq4- (30) 
Ho 

But eqs. (26) and (29) imply /i = 1/2Q, where in this case Q has to be identified 
with the effective mass of the muon pair. Of course, also the first order perturba- 
tive QCD corrections (Compton scattering and gluon production [20]) give rise to an 
asymmetry parameter which, for |qj_| << Q, is again of the type Qj_/Q 2 ; however, 
these perturbative contributions fulfil the Lam- Tung relation[22], which, instead, re- 
sults to be rather strongly violated[20]. Figs. 1 and 2 show comparison between 
formula (30) and the parameter z/, as results from best fits to Drell-Yan data[12]. In 
this connection we observe that, owing result (29), the Q 2 -dependence for DY single 
spin asymmetry with a tranversely polarized nucleon coincides with the one found in 
ref. [24]. 




3 4 5 6 7 8 9 



Q (GeV) 

Figure 2: Azimuthal asymetry in unpolarized DY: the asymmetry parameter v vs the 
effective mass Q of the final lepton pair at fixed |q_i_|. = 23.2 GeV. Data from 
2 nd ref. [20] and fitted by formula (30), taking into account eq. (29). A • = 2.52 
GeV 2 . 

2) CLAS[25] results of SIDIS with longitudinally polarized beam and target. The 
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single spin asymmetry includes a sin2(f) term, characterized by the convolutive prod- 
uct h\ L ® where is the Collins[26] fragmentation function. The function 

hi L (x) = 7r(pl)hi L (x,0) (31) 

has been extracted [25] for li-quarks from that asymmetry, taking into account a model 
evaluation[27] of the Collins fragmentation function and assuming a Gaussian behav- 
ior as regards the p 2 ± dependence of hj~ L (x, p 2 ± ). However the normalization adopted 
in ref. [25] for the "soft" functions is the usual one [2], which, as already explained in 
sect. 2, is different from the one used in the parametrization (6)-(8). In particular, 
as regards h^ L , the conversion factor is C = 2M/Qx, according to eqs. (12) and (29): 
it has to be taken into account, since eqs. (25) hold within our normalization. With 
this re-normalization in mind, we compare (fig. 3) the CLAS[25] results of —xhi^(x) 
[from now on, the "tilde" will be dropped from eq. (31)] to the HERMES data[28] of 
xg±(x), assuming, again, a Gaussian behavior for pj_). The curve corresponds 

to the parametrization of xg\(x) given by GRSV[29] - LO, valence scenario - scaled 
with a factor (1 + i?) _1 [28]. The discrepancy, which appears for x > 0.3, is typical of 
higher twist contribution, not negligible for the modest Q 2 -values involved, 1.5 to 3 
(GeV/c) 2 . 

3) Unpolarized SIDIS coscj) asymmetry, <fi being the azimuthal angle of the final 
hadron momentum with respect to the reaction plane. Recently Anselmino et al.[18] 
fitted EMC [30] and E665[31] data by means of relation (27). They parametrized 
fi(x, p 2 ± ) by making the Gaussian assumption, with (pj_) = 0.25 (GeV/c) 2 , and taking 
the MRST 2001 (LO)[32] parametrization for fi(x). The agreement is good. It is 
to be noticed, however, that the values of x involved are considerably smaller than 
those corresponding to CLAS data, whereas the Q 2 -values are much greater on the 
average. 

6 Conclusion and outlook 

To conclude, relations (24) to (27) have some predictive power, supported by present 
data. This is a goad to further measurements and investigations in forthcoming SIDIS 
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Figure 3: The functions zg" (z)[28] (full circles, Q 2 = 2.5 GeV 2 ) and —xh^[25] (open 
squares, Q 2 = 1.5-3 GeV 2 ). The latter has been re- normalized according to the con- 
version factor C = 2M/Q. The curve corresponds to the GRSV2000 parametrization 
(LO, valence scenario) [29] of xg±(x), scaled with a factor (1 + i?) _1 [28]. 



and DY experiments. 

- First of all, we suggest to determine the sin2(f) SIDIS azimuthal asymmetry 
with a longitudinally polarized target, at higher energies (12 GeV) and at a higher 
precision than in ref. [25]. In particular, new extractions of hj^ L (x) should take into 
account the recent determination of the Collins function by BELLE [33] collaboration. 
According to our result, this asymmetry is predicted to decrease like Q~ 2 . Moreover 
we could test our prediction about higher twist effects, according to which the sum 
hf L (x) + gi(x) decreases again like Q~ 2 . The function h^ L (x) could be determined 
also by means of the DY single spin asymmetry, which includes a term of the type 
h\L ® h\\ in this case it would be suitable to realize pp collisions [34, 35] or n~p 
scattering [36], since the function hf can be extracted from unpolarized DY. 

- The above determinations, while interesting in themselves, would increase relia- 
bility of the approximate relations deduced above, in particular, of the last two eqs. 
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(25), not directly testable at the moment, which can be exploited for determining 
approximately hy? and h^ T through the chiral-even function g 1T . Indeed, this may 
be extracted from double spin asymmetry in SIDIS, with a longitudinally polarized 
beam and a transversely polarized target [37]. No particular problems should arise in 
inferring gyr from asymmetry data, since this asymmetry is characterized by the con- 
volutive product giT<S>D, where D is the common unpolarized fragmentation function 
of the pion. 
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